Approximations of permutation-symmetric vertex couplings 

in quantum graphs 

Pavel Exner and Ondfej Turek 

Abstract. We consider boundary conditions at the vertex of a star graph 
which make Schrodinger operators on the graph self-adjoint, in particular, the 
two-parameter family of such conditions invariant with respect to permuta- 
tions of graph edges. It is proved that the corresponding operators can be 
approximated in the norm-resolvent sense by elements of another Schrodinger 
operator family on the same graph in which the <5 coupling is imposed at the 
vertex and an additional point interaction is placed at each edge provided the 
coupling parameters are properly chosen. 



1. Introduction 

There is no necessity to describe here in extenso what quantum graphs are 
and why they are important; if such a need nevertheless arises we can refer to 
papers from the daw ns of th e history RS53 , from the times of new beginning in 
the eighties [GP88L IES8 9: . to more recent work containing a rich bibliography 
KS9Q1 IKu04j . and last not least, to the other contributions making this volume. 

As in the most of the mentioned work, the object of our interest here are 
Schrodinger operators on metric graphs; we neglect external fields and consider a 
free spinless particle on the graph, with the Hamiltonian which acts as Htpj = —ip'j, 
where ipj denotes the wave function at the jth edge. It is known for longtime 
ES89 that in order to make H self-adjoint, a vertex joining n graph edges may be 
characterized by boundary conditions involving n 2 real parameters; they have the 
form of a linear relation between W(0), the column vector of the boundary values at 
the vertex (identified conventionally with the origin of the coordinates), and ^'(O), 
the vector of the derivatives, taken all in the outgoing direction. 

A general and elegant form of these boundary conditions was found in KS99 : 
any vertex in which n edges meet can be described by a pair of n x n matrices A, B 
such that rank (A, B) = n and the product AB* is self-adjoint. The self-adjointness 
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is guaranteed if the corresponding boundary values satisfy the condition 

(1.1) A*(0) + £#'(0) = . 

Moreover, soon after several authors |FTOOL IHaOOL IKSOOj pointed out that the 
matrix pair in l|l.l|) can be made unique by choosing 

(1.2) A=U-I, B = i{U + I) 

with a unitary U ; a nontrivial coupling between the edges corresponds naturally to 
the situation when the matrix U is non-diagonal. A simple proof of this fact for 
n = 2 was given in |FTOO| and extended to any n in |CE04| . 

While the conditions (|1.1|) ensure self-adjointness of quantum graph Hamiltoni- 
ans, or in physical terms conservation of probability current in the vertex, they say 
nothing about a physical meaning of such a vertex coupling. A natural way to ad- 
dress the last question is to investigate approximations of a quantum graph by more 
realistic systems with no free parameters. An example is a quantum particle living 
in a configuration space in the form of a thin tube-like domain; one can consider a 
family of such domains shrinking to the given graph. A solution to this problem at 
the level of eigenvalue convergence was found [KZOll IRSOll ISaOlL IEP05j in the 
situation that the tube-like domain supports Laplacian with Neumann boundary 
conditions (or similar operators), and an extension to the resolvent convergence 
has been announced |Po05j . These results, however, gave a partial answer to the 
problem stated above because the limit leads to the free boundary conditions, 

n 

(1.3) ^(0) = Vfc(0), j,k = l,...,n, ^>;-(0) = 0, 

j=i 

only. It is hoped that other approximating families, say, using Dirichlet Lapla- 
cians, geometrically induced and/or external potentials, could yield different vertex 
couplings, but this problem is difficult and no such results are known at present. 

A less ambitious program aims at approximating vertex couplings by means of 
Schrodinger operators on the graph itself, using suitable families of scaled potential, 
regular or singular. This is relatively easy as long as we attempt to approximate 
couplings with wavefunctions continuous at the vertex, i.e. the one-parameter fam- 
ily of the so-called 6 couplings Ex96a 4 described by the conditions 

n 

(1.4) ^(0)=V*(0)=:iK0), j,k = l,...,n, £ $(0) = c^(0) 

3=1 

with a S M, which is obtained from Ijl.lfl and (|1.2|) by choosing U = n + ia J — I, 
where J and I are the n x n matrix whose all entries are equal to one and the 
unit matrix, respectively. The procedure is analogous to the approximation of 5 
interaction of the line AGHH : one starts with the conditions l|1.3fl and adds at 
each edge a naturally scaled potential, Vj >e (x) = -Vj (-J for some Vj S L 1 . A 
norm-resolvent limit then yields the S coupling with a := • J Vj(x) dx [Ex96b . 

The situation is more complicated if the wavefunctions are discontinuous at the 
vertex. The simplest example of such a situation is the so-called 5^ coupling, 

n 

(1.5) $(0) = V£(0)=:^(0), i,fc=l,...,n, £>i(0) = /ty'(0) 

3=1 
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with (3 £ R. An inspiration can be found in the way in which Cheon and Shigehara 
CS98a, CS98b approximated formally the 5' interaction on the line using a non- 
linearly scaled family of S interactions - their argument was later shown to yield 
a norm-resolvent convergence and to lead to approximations in terms of regular 
potentials [ANOOl IENZ01| . It was shown in CE04 1 that the CS-type method 
can be used to approximate the 5' a coupling for any n, the approximating oper- 
ator domains having functions continuous at the vertex. The aim of this paper 
is to show that this result can be extended to couplings with discontinuous wave- 
functions which are invariant with respect to permutations of the graphs edges: an 
approximation using a S coupling at the vertex and an n-tuple of 5 interactions at 
the edges approaching the vertex will be derived for all such couplings 2 . We will 
see that in the generic case the idea of CS98a has a direct, albeit rather tedious 
extension to the graph case, while for two one-parameter subfamilies the choice of 
coupling parameters requires a modification. Extensions to more general boundary 
conditions inspired by SMC99 and approximations by regular potentials are left 
to a subsequent publication. 

2. Permutation symmetric vertex couplings 

As in the previous work cited above we consider a star graph T consisting of 
n halflines meeting at a single vertex. The corresponding Hilbert space is thus 
H = L 2 (R + ). A general Hamiltonian describing a free particle living on 

the graph is a self-adjoint extension of the operator Hq acting as Hgipj — —ip" on 
functions * = {ipj} G 0™ =1 W 2 - 2 (R + ) satisfying the conditions ^-(0) = ^.(0) = 0; 
each such extension is specified by a boundary condition (| 1 . If) at the vertex. 

Since the action of these operators at each component of the wavefunction is the 
same, symmetry properties of the extensions are given by those of the boundary 
conditions. We will be interested in the permutation-invariant extensions, first 
introduced in |ES89 , which form a two-parameter family. 



PROPOSITION 2.1. The boundary conditions ^1.1)) are permutation invariant if 
and only if the matrix U in equals 

(2.1) U = aI + bJ 
with complex coefficients a, b satisfying the relations 

(2.2) \a\ — 1 and \a + nb\ = 1. 

Proof. The condition (|l.lfl is permutation invariant iff it is satisfied at the 
same time by the vectors P\I/(0) and PvE^O) for any P G S n . Multiplying it by 
P -1 from the left we get 

{P- 1 UP - J)*(0) + i{p- x VP + 7)*'(0) = , 

the matrix P~ 1 UP being obviously unitary. In view of the uniqueness of the 
parametrization (|1.2fl the property is equivalent to P~ 1 UP = U. Next we notice 
that a simultaneous permutation of the rows and columns leaves the diagonal ele- 
ments on the diagonal, and the off-diagonal ones off the diagonal; since P~ X UP = U 



small correction is due with respect to this paper: in the statement and proof of Theorem 1 
/3/n should be replaced everywhere by j3. 

2 We will be concerned with nontrivial couplings only corresponding to non-diagonal matrices 
U; in the opposite case the task is reduced to a much simpler halfiine problem - cf. fFCT02 . 
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has to be satisfied for any P £ S n it follows that U = al + bJ for some a, b G C. 
The conditions restricting the values of a,b follow from the unitarity of U, 

{UU*) l3 = \a\ 2 5 l3 + 25i( a 5) + n\b\ 2 = % , 

which yields the relations 

H 2 + 23?(a6) + n|6| 2 = 1 , 25R(a6) + n\b\ 2 = , 

for i = j and i =/= j, respectively. Substituting from the second to the first one we 
get \a\ 2 = 1. Finally, using \a + nb\ 2 = \a\ 2 + 2n$i(ab) + n 2 \b\ 2 we see that the 
left-hand side of the second relation is a multiple of |a + nfo| 2 — |a| 2 . □ 

For dcfinitcncss we will denote in the following the self-adjoint extension cor- 
responding to fixed a, b as H a b . Notice that the boundary conditions described by 
Proposition 12 . II can be also written more explicitly as the following system, 

n n 

(a - 1)^(0) + bJ2 MO) + i(a + l)^-(O) + i&X>U0) = , j = l,...,n, 

k=l k=l 

which shows, in particular, that (a — X)ipj{0) + i(a + 1)^(0) is independent of j. 
To get a useful equivalent formulation we subtract the fcth one of these condition 
from the jth one obtaining 

(2.3) (a - 1) (V^-(O) - ^ fc (0)) +i(a+ 1) ($(0) ~ V4(0)) = , j, k = 1, . . . , n , 
while summing all of them gives 

n n 

(2.4) {a-l + nb)J2 M°) + «(«+! + nb) ^ ^(0) = . 

fc=i fc=i 

Examples 2.2. We have already mentioned that a = — 1 and b = 2 . de- 
scribes the 5 coupling (|1.4|l . similarly a = 1 and b = i ^'_ n corresponds to the S' s 
coupling HI. 51) . Another example is the 8' coupling Ex96a , 

(2.5) ^^(0) = 0, ^(0)-^(Q) = ^(0)-^fc(0)), j,k = l,...,n, 

referring to a = and b = n 2 ^ , and its dual counterpart S p with the roles of 
functions and derivatives interchanged for which a = "7 M and b — jv- . 

3. Approximation: a heuristic argument 

Let us describe the family we will employ to approximate permutation-symmet- 
ric Hamiltonians H a,b . Let us recall that we will consider all operators of this class 
with the exception of those with a 6 coupling, i.e. with the wavefunctions continuous 
at the vertex, because for the latter we have the natural approximation described 
in the introduction. We denote by H UtV (d) the operator which is obtained from 
H u ^ := /f- 1 ) 2 /(™+ l «) by adding a 8 interaction of strength v to each edge at the 
distance d from the vertex; it is the same scheme which was used in the particular 
case of 5' s treated in CE04 . The aim of the present section is to derive formally 
how the values of the parameters u, v as functions of d should be chosen. 

At the vertex the boundary condition defining H U}V (d) are of the form (|1.4|) 
with a replaced by u, while the added 5 interactions are characterized by 

(3.1) ^(d+)=Vj(d-)=:^(d), i>' j (d+)-TP' j (d-)=vip j (d), j = l,...,n. 
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To find relations between the boundary values, we employ Taylor expansion 

(3.2) ^(d)=^(0)+d^(0) + O(d 2 ), ^(d-)=^(O+) + 0(d), j = l,...,n; 

we want to choose u, v to get the relations Ij2.3|l and (|2.4(l in the limit d — > 0+. The 
first one of the relations (|3.2|1 together with the continuity at the vertex imply 

(3.3) ^(d) - Md) = d($(0) - Vfe(O)) + 0{d 2 ) . 

Furthermore, the second one of the relations (|3.2|) in combination with 13.1fl tell us 
that the difference ipj(0+) — ip' k (0+) is equal to 

$(d-) - i>' k (d-) + o(d) = ^(d+) - ^(d+) - «(^(d) - ^ fc (d)) + o(d) 

giving thus d (ipj(d+) - i/)' k (d+)) - v(i/)j(d) - ipk{d))) + 0(d 2 ) as the value of the 
left-hand side in (|3.3I) . which can be rewritten as 

(f + dv) (Vv(d) - Md)) - d (^(d+) - ^(d+)) = 0(d 2 ) • 

This should give (a - 1)(^(0) - Vfe(O)) + i(a + l)(ip' j (0+) - tp' k {°+)) = in the 
limit d — > 0+. As we have mentioned above, the case of a <5 interaction in which 
we have a = — 1 is excluded, hence we are allowed to require 1+ J* V — ^"'^ . This 
in turn yields the following relation for the parameter v, 

, x 1 a - 1 

3.4 V = — -{ -; 

y ' d a + 1 

notice that it is real- valued in view of the condition \a\ — 1, because 

a- I \a\ 2 + 2i^a-l 3a 

i = i^—, m = -2- ^ 61. 

o+l |a+l| 2 |a+l| 2 

It remains to find u. We employ again the first of the relations H3.2J) together with 
both the vertex conditions (|1 .4|> for a = u rewriting in this way —1 ipj(d) as 

n n n 

n 



W(0) + dJ2^j(0+) + 0(d 2 ) = -^^(O+) + dJ^(O+) + 0(d 2 
/ 1 " i=i i=i 

As before we use l|3.1f) and l|3.2|l to eliminate ^-(Q+), 



E $ (0+) = E $ (<*-) + o(d) = J2 $ - w E ^ ( d ) + °( d ) 

3=1 j =1 i =1 i =1 

Substituting into the expression for ipj{d) we get after a simple manipulation 



(1 +« £ + d)) E^w = £ +rf) E + w + 

using the value of v given by l|3.4|l we find that the quantity 



^tt) ; + 'fir-) £ *« + (;+«) (t + ow) 



6 



PAVEL EXNER. AND ONDREJ TUREK 



behaves as 0(d 2 ) in the limit d — ► 0+. We will look for u having a stronger 
singularity than v assuming ^ = 0(d 2 ); then the last claim simplifies as follows, 



In .a — 1 



d 2 



v 



l ^rj j + E#+) = 0(d) . 



3=1 3=1 

This is required to give the condition (|2.4|) in the limit d — > 0+ which happens if 

In a— 1 a— l + 



<i 2 it a + 1 i(a + 1 + n&) ' 

provided the two denominators containing the coupling parameters do not vanish. 
The first one is zero for the 6 coupling which we have excluded from the outset, the 
second one vanishes iff a, b correspond to the 8 P coupling described in Examplcs l2.2l 
It is also clear that in view of the conditions \a\ = 1 a |a + n&| = 1 the fractions 
°+i+"b an d aTT are P ure ly imaginary. This motivates us to choose 

. n fa - 1 + nb a - l x 

(3.6) u = i— H - 

y ' d 2 \a + 1 + nb a + 1 

assuming that the expression in the parentheses is nonzero which is true as long as 



(3.7) a(a + nb)^l. 

The parameter u defined by Q3.6JI is, of course, real and u = 0(d~ 2 ) as d — > 0+; 
this concludes our search for the approximating operator family in the generic case. 

It remains to carry on the heuristic argument for the two excluded one-parame- 
ter subfamilies, the S p coupling and the one violating the condition (|3.7|l . We will 
show that the coupling of the 6 interactions at the graph arms can be preserved, 
it is only necessary to change the function u describing the vertex. Let us first 
suppose that the latter has a stronger singularity at d = 0, for instance, 

(3.8) u=^ 

for a fixed nonzero £ S R (in fact, one can replace d 3 by d v for any v > 2). 
Substituting this into i|3.5|) we get a condition which in the limit d — > 0+ yields 



a 
i — 



1 n n 



The left-hand side makes sense since a ^ — 1 and it is easy to check that if (|3.7|l is 
not valid, i.e. a + nb = a -1 , the last relation is equivalent to l|2.4|l . On the other 
hand, to deal with the S p coupling we take u with a pole singularity, 

(3.9) u = ~. 

The second term in (|3.5|) then vanishes and we find that 53j=i V'i(O) = 0(d 2 ) which 
gives in the limit d — » 0+ the condition l|2.4|l for the particular case of 8 p . 
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4. The main result 

Now we are ready to formulate and prove our main result. 

Theorem 4.1. Given complex numbers a —I and b satisfying the con- 
ditions 12.2(1 . define u — u(d) and v = v(d) for d > as in the previous section, 
i.e. by the relations ([3.6(1 . ((3.8(1 . (13.9(1 . and (13.4(1 . respectively; then the operators 
H UtV (d) converge to H a,b in the norm resolvent topology as d — > 0+. 

PROOF. To begin with we observe that the permutation symmetry of the 
boundary conditions 1(2.3(1 and 1(2.4(1 allows us to simplify the task by reducing 
it to independent halfline problems. To this aim let us find the spectrum of the 
matrix U . Since the latter equals U = al + bJ it is sufficient to look at the matrix J 
which has rank one, and thus zero is its eigenvalue of multiplicity n — 1 correspond- 
ing to vectors with vanishing component sum; the remaining simple eigenvalue is 
n. The corresponding eigenvalues of U are a and a + nb, respectively, with the 
same multiplicities. Consider first the generic case where these eigenvalues are not 
inverse to each other by 1(3.7(1 and none of them equals to —1. 

The symmetry allows us to decompose the operator in question, H a ' b on the 
Hilbert space H = 0™ L 2 (R+) = L 2 (R+) ® C", into orthogonal sum of two com- 
ponents. The first one denoted as H^ a ' b acts at the "scalar" subspace isomorphic 
to i 2 (M + ) ® C consisting of functions if? G H which are symmetric with respect 
to permutations, ^j{x) — ipk(x) for all j,k = 1, ...,n; it is characterized by the 
boundary conditions (a + nb — 1)^(0) + i(a + nb + 1)^'(0) = 0. The other one for 
which we use the symbol fj( n - 1 ) a > b acts on the orthogonal complement which is 
isomorphic to L 2 (R + ) ® C™ -1 consisting of \& G Ti with vanishing component sum. 
The action of Jj( n - 1 ) a < b on all linear combinations Y^j=i c jV'j( a; ) is identical and 
the boundary conditions arc (a - 1)^(0) + i(a + 1)*'(0) = 0. 

In the same way one can decompose the approximating operators H UtV (d). 
The part Hu]l(d) acts on the "scalar" subspace of functions invariant with respect 
to permutations, the boundary conditions being ^'(0) = ^^(0). The remaining 
component acts on "(n-l)-dimensional vector functions" being isomorphic to n — 1 
copies of the "scalar" problem with Dirichlet boundary conditions. 

We will use the fact that the resolvents of all the involved operators can be 
constructed explicitly using a standard ODE result in combination with Krein's 
formula. Let us consider first the part independent of the coupling at the vertex. 
For a fixed k from the upper complex halfplane the Green function of the Laplacian 
on the halfline with Dirichlet condition at the origin is 

Gk(x,y) = ysin(kx < ) ikx > = - sinh(ra<)- Ka;> , 
k k 

where we denote conventionally x < — min{x,?/}, x > = max{x,y}, and n = —ik. 
The S interaction at the point x = d represents a rank-one perturbation of the 
above free resolvent, and corresponding Green's function is found easily with the 
help of (EH) as in AGHH Sec. 1.3] or CE04] to be equal to 

,s r v ( s r i s , Qk(x,d)g k (d,y) 

4- 1 Gk(x,y) = Gk(x,y) + — — — ■ 

-v 1 - y k {d,d) 

Next we have to find the Green function of the approximated operator. Following 
|Wel Sec. 8.4] we need a solution of the equation — rjj" = k 2 ip satisfying the condition 
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(a - l)V'(O) + i(a + l)ip'(0) — and its Wronskian with <f>(x) — e KX ; this yields 

(i(a — 1) sinh/ta^ + n(a + 1) coshKx<) e~ KX> 
(42) G ^ V) = «(i(a-l) + «(a + l)) ' 

Now we are going to show that Q\ K converges to Q? K pointwise as d — > 0+. First we 
suppose that both arguments are not smaller than d; without loss of generality we 
may put d < x < y rewriting the difference G\ K {x, y) — Gf K (x, y) as 

~ Ky I , e~ KX sinh 2 Kd i(a— 1) sinhfta; + K,(a+1) cosIikx 

sinh kx + 



^j 1 ' -e-^sinh/td i{a-l) + K(a+1) 

a+l+id(a—l) 

The sum of the first and the third term in the bracket equals i(a-i)+l(a+i) being 
independent of d. In the second term we use the expansion sinh(a;) = x + 0(x 3 ) 
and obtain after simple manipulations 



— KX 



sinh 2 Kd -kx f «(a+l) 



^j 1 ' -e-^sinh^d \i(a-l) + K{a+l) 

a-\-l-\-'ia(a— 1) 



+ 0(d 2 ) 



hence the terms non- vanishing in the limit cancel and there is a K > such that 

\gv K (x,y)-gs t (x,y)\<Ke- K *e-"y d 2 

holds for all d < x < y, and by the same argument also for d < y < x. 

Next we suppose that x < d < y when the Green's function difference is 

~ Ky ( . , sinh Kxe~ Kd sinh 2 Kd i(a — 1) sinh^a; + n(a+l) cosIikx 

1 sinh kx ■ 1 



-e-^sinh/td i(a-l) + K{a+l) 

We want to show that the expression in the bracket is uniformly bounded in x, d 
provided d is small enough; we may suppose that d < 1. The first and the third 
term are bounded in view of the continuity at x = 0, the middle one is easily found 
to be 0(1) as d — > 0+. Hence there is an L > independent of x,y and d such 
that \Gi K (x,y) — Gi K {x,y)\ < Le"^^ 9 , and in a similar way one can estimate the 
Green function difference for y < d < x. It remains to deal with the case when 
both argument are less than d, say x < y < d; we may again suppose that d < 1. 
The middle term in the above expression is then replaced by 

sinh kx sinh Ky e~ 2Kd 

Kd{a+l) _ Kd • , j 
a+l+id(a-l) C Sml1 KCl 

and one checks easily that there is an M > independent of x, y and d such that 
the pointwise estimate \Gf K (x,y) — Gf K {x,y)\ < M holds. 

These bounds allow us to estimate Hilbert-Schmidt norm of the difference, 



+ 



-- / / \G? K (x,y)-G? K (x,y)\ 2 dxdy 
Jo Jo 

pd poo 

GUx,y)~G? K (x,y)\ 2 dxdy+ / \Gl(x,y)-G? K (x,y)\ 2 dxdy 

Jo Jd 

poo pd poc poo 

/ / \GUx,y)~G: K (x,y)\ 2 dxdy+ / / \gV K {x,y)-g? K {x,y)\ 2 dxdv 

Jd JO Jd Jd 
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It is straightforward to check that last integral does not exceed the value 

f>OC/>00. .0 / ts- \ 2 



and similarly the first one and the middle two are estimated by 4M 2 rf 2 and 2 $t( R ) d, 
respectively, which means that 

2 

= 0. 



lim 



2 



2 



and the same is a fortiori true for the operator norm. This concludes the argument 
for the first component of the operator. 

The proof for the "scalar" component is similar, just a bit more complicated, so 
we can skip some details. First we construct the Green function for the 6 coupling 
with the parameter agl projected on the subspace of functions with coinciding 
components; in a similar way as above we find that the resolvent kernel equals 



e > (u 
-Ax,y) = —7^— — r ( - £ 



sinh kx< + K cosh / 



An analogous construction for the (negative) Laplacian with the boundary condition 
(a + nb — l)V'(O) + i(a + nb + l)i//(0) = at the origin gives the Green function 

a .b, e~~ Ky (i(a + nb — 1) sinh kx + re(a + nb + 1) cosh kx) 

lK ^ X,V > = K (i(a + nb~- 1) + K {a + nb+ 1)) ' 

Finally the resolvent kernel of the approximating function, to be compared with 
the last expression, is obtained again from Qf K (x,y) by means of Krein's formula 

r u,v, , ru i x , d)Qf K (d,y) 
To estimate the Green function difference we assume again first that d < x < y, so 

/~>U,V / \ s-rCL,b / \ 

Sin fay) - Gin ( x >y) 

,2 



„„ /„ . , , ~r~ — 77 (~ sinh Kd + k cosh nd) 

e~ Ky / - smhnx + Kcoshra Vn 1 



« \ 1 + k ^iS.-d - IS (3 sinh Kd + K cosh Kd ) 

i(a + nb — 1) sinh kx + k{cl + nb + 1) cosh kx \ 
i(a + nb - 1) + «(a + nb + 1) J 

where we have used (13. 411 and u should be substituted from l|3.6|l . Our aim is to 
find the behavior of this expression for small d using the expansion 

cosh(x) = l + 0(x 2 ), — !— = 1 - x + 0(x 2 ) asx->0. 
1 + x 

The first term gives 

— sinh kx + k cosh kx „ 
* = sinh kx + e~ KX 0(d 2 ) , 

n + K 

for the second one we get after a straightforward but tedious computation 

k{cl + 1 + nb) 



k(ci + 1 + nb) + i(a — 1 + nb) 



0(d) 
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and the third is independent of d; putting everything together we find 

G^(x,y)-g:*(x,y)=e-™e- K yO(d), 

because the non-vanishing terms cancel again. In other words, there is a K' > 
independent of x, y and d such that the following inequality 



(4.3) 



< K'e- KX e~ Ky d 



holds. In the "mixed" < d < y, we have 



Q l k (z> v) - {x, y) = 



■ sinh kx + k cosh kx 



- + 



+ 



-^-j (- sinh kx + k cosh kx) (- sinh Kd + k cosh Kd) 



(*+*) 



d(a+l) 



r-p- ( - sinh Kd + k cosh Kd) 

-+k \ n ) 



n a+\+id{a-l) 

i(a + nb — 1) sinh kx + k(cl + nb + 1) cosh kx 
i(a + nb — 1) + k(cl + nb + 1) 

The first and the third term at the right-hand side are obviously bounded indepen- 
dently of x, y and d, and in the same way as above one can check that the second 
one is 0(1) as d — > 0+, hence there is an L' independent of x,y and d < 1 such 
that 



s->u,v I \ /->a.b r \ 



The same is naturally true if the roles of x and y are interchanged. It remains to 
analyze the situation when both x, y do not exceed d, say x < y < d, when 



e Ky % sinh kx + k cosh kx 



+ 



±+K 

( - sinh kx + k cosh kx) ( - sinh ny + k cosh Ky) 



rrfer^ n — %t— f — sinh nd + k cosh Kd) 

a+l+id(a— 1) ^ n ' 

e _KS ' i(a + nfe — 1) sinh kx + k(o, + nb + 1) cosh kx 
k i(a + nb — 1) + k(cl + nb + 1) 

In the same way as above one establishes existence of an M' > independent of 
x, y and d < 1 such that 



S^(x,y)-g^(x,y) 



< M' . 



Using these bounds and repeating the above Hilbert-Schmidt estimate we get 



lim 

d^0+ 



^H { Jl{d)^ 2 } ~ ^H(i)«.l> (^ 2 ) 



0. 



which implies the analogous limiting relation for the operator norm of the resolvent 
difference which we set out to prove. 

In the remaining two cases it is sufficient to consider the "scalar" component 
because the orthogonal complement does not contain the parameter u. Take first 
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the case when the condition (|3.7|l is violated. If the variables satisfy d < x < y we 
can rewrite the Green function difference using (|3.8(l as 



expanding the first two terms at the right-hand side we establish existence of a 
K' > independent of x, y and d such that the inequality (|4.3(l holds. In a similar 
way one proceeds when one or both arguments are smaller than d. The same can 
be done in the S p case where the resolvent difference for d < x < y is 



and the other variable combinations are dealt with analogously. The Hilbert- 
Schmidt estimate is the same as in the generic case; this concludes the proof. □ 



We have mentioned in the introduction that approximation including singular 
couplings can be used an intermediate step in a search for approximations based on 
regular potentials. In this sense S coupling and S interactions are preferable because 
in this case we already know how to make the second step; hence our result paves 
way to a complete potential approximation of permutation symmetric couplings. 

In particular, comparing with jCE04| we do not need S p coupling to approx- 
imate 5' , and the S p itself can be approximated by S interactions. We have seen, 
however, that in this case the central singularity is of a pole type with respect to 
d similarly as the couplings of the 5's at graph edges. This illustrates the excep- 
tional character of 5 P which is in a sense akin to 5, with the roles of the "scalar" 
and (n — l)-components interchanged. The remaining one-parameter family of cou- 
plings violating the condition (|3.7|) needs, on the contrary, a stronger singularity 
with respect to d at the vertex. The reason of this behavior is not clear; this un- 
derlines one more time the fact that our present understanding to the zoology of 
vertex couplings in quantum graphs is still far from satisfactory. 
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